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We present a comprehensive theoretical study of the resonance fluorescence spectra of an exciton- 
driven quantum dot (QD) placed inside a high-Q semiconductor cavity and interacting with an 
acoustic phonon bath. We derive a quantum master equation (ME) in the polaron frame which 
includes exciton-phonon and exciton-cavity coupling to all orders. This work details and extends the 
theory used in a recent issue of Physical Review Letters (C. Roy and S. Hughes 20f 1: Phys. Rev. Lett. 
106 247403) to describe the QD Mollow triplet in the regime of semiconductor cavity-QED. Here we 
introduce two ME forms, Nakajima-Zwanzig (NZ) and time-convolutionless (TC), both to second 
order in the system-phonon-reservoir perturbation. In the polaron frame, these two ME forms 
are shown to yield equivalent population dynamics and fluorescence spectra for a continuous wave 
(cw) driving field. We also demonstrate that a Markov approximation is vaUd for computing the 
incoherent scattering processes and we subsequently exploit the Markovian TC ME to explore the 
resonance fluorescence spectra of an exciton-driven QD. Both cavity-emitted and exciton-emitted 
spectra are studied and these are found to have qualitatively different spectral features. Using a 
coherent driving field, the well known characteristics of the atomic Mollow triplet are shown to 
be considerably modified with electron-acoustic-phonon scattering and we highlight the key effects 
arising from both cavity coupling and electron-phonon coupling. Regimes of pronounced cavity 
feeding and anharmonic cavity-QED are exemplified, and we find that the cavity coupling depends 
sensitively on the exciton-cavity detuning and the temperature of the phonon bath. We show how 
the full width at half maximum (linewidth) of the Mollow triplet sidebands varies as a function of 
the square of the Rabi frequency of the cw pump. Phonon-mediated cavity coupling also contributes 
to the spectral broadening of the Mollow triplet, depending upon the exciton-cavity detuning and 
the strength of the exciton-cavity coupling rate. Finally, we calculate the fluorescence spectra for 
off-resonance cw driving and investigate the resulting MoUow-triplet linewidths. 

PACS numbers: 42.50.Ct, 78.67.Hc, 78.55.-m 



I. INTRODUCTION 

Semiconductor quantum dots (QDs) interacting with 
microcavities are a field of intense research where 
ideas from conventional cavity-quantum electrodynam- 
ics (cavity-QED) are being applied and studiecP'?. 
These chip-based semiconductor devices show enormous 
promise for technological applications ranging from quan- 
tum information processing and quantum metrology to 
the creation of single photon on demand^ ^. The pos- 
sibility of exploiting quantum anharmonicities in these 
cavity-QED systems may also lead to photon blockade^^. 
Additionally, the spin degree of freedom of the underly- 
ing exciton states is of interest for quantum information 
processings^. In the past few years, an increasing number 
of experimental studies have focussed on resonance fluo- 
rescence of a QD coupled to a cavity mod^^^HUl, Progress 
has been made in the study of off-resonant QD cavity sys- 
tem that can be used to observe the photoluminescence 
(PL) spectra of the QD or the cavitjffsffisl. 

Recent experimental studies of resonance fluorescence 
of InGaAs QDs embedded in a high-quality micropillar 
cavity^ have provided clear demonstration of excitation- 
induced dephasing (EID). In QD structures with no 
cavity coupling, the phonon-induced EID process is re- 
sponsible for the intensity damping of the QD exci- 
tonic Rabi rotations induced by pulsed excitatiorPsMl, 
The experimental resonance fluorescence spectra in QD- 



cavities is similar to the characteristic Mollow triplet of 
an atomic system, but the sidebands of the triplet are 
found to undergo systematic spectral sideband broaden- 
ing with increasing continuous-wave (cw) drive^l. This 
drive-dependent broadening results from the interaction 
of the driven QD with the underlying phonon rcscr- 
For cavity-QED (quantum electrodynamics) 
systems, acoustic phonon scattering processes also intro- 
duce an additional co upling mechanism between a cavity 
mode and the excitorPSEH. This phonon-mediated cav- 
ity coupling depends on the QD-cavity detuning and the 
intrinsic parameters of the phonon bath. Cavity-exciton 
coupling can be used to measure the fluorescence spec- 
tra of the QD via the cavity mode which is convenient 
because of spectral (and possible geometrical) separation 
between the exciton and the cavity mode. 

Acoustic phonon processes also play an important role 
in incoherently pumped cavity-QED systems which are 
known to result in off-resonant "cavity feeding"'^2H22l^ ^n 
effect which also exists ( and i s particularly strong) in 
coherently driven systemd28129|_ There have been sev- 
eral attempts at an improved theoretical description of a 
semiconductor cavity-QED system to include the effects 
of electron-acosutic-phonon scatterin^^IISBIlII ^j-^ 

ample, the vacuum Rabi doublet of a cavity-QED sys- 
tem in the strong coupling regime undergoes modification 
due to phonon coupling^, e.g., the spectral doublets — 
for on-resonance coupling — are asymmetric, especially at 
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low temperatures, which can be e xplained in terms of 
phonon-assisted processe^^^'^^'^^EIl, 

The complexities of modelling a coherently-driven 
semiconductor cavity-QED system (where multiphoton 
effects are required) are best treated in a quantum mas- 
ter equation (ME) formalism which can be used to inves- 
tigate both qualitative and quantitative changes in the 
Mollow spectra due to phonon-induced scattering. In 
this paper, we present the theoretical details behind the 
model used in a recent edition of Physical Review Let- 
ters^^. We then use this model to systematically explore 
the exciton-driven Mollow triplet in a wide range of ex- 
citation regimes. We first present two alternative ME 
forms which govern the evolution of the reduced density 
operator, wherein we incl ude the s ystem-bath incoherent 
interaction to second ordei'^^'^SMi system-reservoir 
coupling. These two MEs are termed the Nakajima- 
Zwanzig (NZ) form and the time-convolutionless (TC) 
form. To second-order in the polaron transformed 
system-phonon-reservoir coupling, these ME forms are 
equivalent to the standard Born and Born-Markov MEs, 
but we keep the former notation (TC) to avoid confusion 
with a second-Markov approximation which can be ap- 
plied on the memory relaxation of the phonon reservoir. 
The polaron frame, obtained following a unitary polaron- 
transformation, allows us to study the effects of phonon 
dephasing on the coherent part of the Hamiltonian to all 
orders, i.e., beyond the Born approximation. Moreover, 
these polaron MEs allow one to study the nonperturba- 
tive dynamics of a cavity-QED system interacting with 
a radiation reservoir and driven by an external cw laser 
field which is also under the influence of the phonon en- 
vironment. The use of the polaron frame eliminates the 
QD-phonon coupling in the system Hamiltonian and in- 
troduces a modified exciton-cavity coupling and a modi- 
fied radiative decay rate'*^ ''^ . 

The rest of our paper is organized as follows. In Sec.[TT| 
we present the model Hamiltonian and derive a second- 
order TC form of the ME and provide a comparison with 
the second-order NZ form. The computational effort nec- 
essary to solve the NZ ME is significantly more compli- 
cated than the TC ME. However, for cw drives, we show 
that both ME forms produce the same system dynamics 
and spectra. We also study the Markov limit of the TC 
ME by neglecting the memory dynamics of the phonon 
reservoir, and show that, for our parameters and drives, 
the full non-Markovian solution is identical (or at least in- 
distinguishable) to the dynamics obtained in the Markov 
limit. We subsequently use the Markov form of the TC 
ME form for our remaining calculations. In Sec. |III| wc 
present and discuss our numerical results of the resonance 
fluorescence spectra of a coherently driven cavity-QED 
system with an exciton pump field. The QD linear ab- 
sorption spectrum is first shown for several different bath 
temperatures. We then study the role of the temperature 
of the underlying phonon reservoir on the fluorescence 
spectrum. We explore the effects of QD-cavity detun- 
ing and identify regions where the QD-cavity coupling 



is particularly strong. Next, we provide an analysis of 
the effects of the cw-laser drive and QD-cavity coupling 
on the spectrum and present detailed simulations of the 
fluorescence spectra for various QD-cw laser detuning. 
We subsequently investigate EID effects by numerically 
extracting the spectral linewidths of the Mollow triplet 
sidebands for various cw-laser drives. In addition, we 
study the dependence of QD-laser detuning on the Mol- 
low triplet broadening. In Sec. IV 



we summarise our 



results and present our conclusions. 



II. THEORY 

A. Model Hamiltonian 

We start by introducing a model Hamiltonian which 
describes the exciton-driven cavity-QED system, where 
the QD interacts with a phonon reservoir and the dot is 
driven by a classical cw-laser field. Considering a single 
exciton (electron- hole pair), and transforming to a frame 
rotating with respect to the laser pump frequency, Wi, 
the model Hamiltonian is 



H 



where the phonon bath is represented by a continuum 
of harmonic oscillators of frequency with annihila- 
tion (creation) operators bq{b^^), and Xq are the exciton- 
phonon couplings (assumed real); also introduced in 
Eq. ([l]) are the detunings of the exciton (ujx) and cav- 
ity (iOc) from the pump laser frequency, A^l = '^a ^ 
[a = x, c), the cavity mode annihilation (creation) opera- 
tor, a{a) ), the Pauli operators of the exciton, cr"*", and 
the exciton-cavity coupling rate, g. The exciton pump- 
ing rate, r|x^ describes the coherent exciton pumping from 
the cw laser field. We have included the role of phonons 
at the level of the independent boson model (IBM) 
which is known to yield good agreement with experi- 
ments on solid-state QDpl. A more complete descrip- 
tion of the phonon reservoir, e.g., to allow broadening 
of the zero phonon line (ZPL), would require the gen- 
eralization of the above system Hamiltonian to include 
additional microscopic pr ocesse s, such as spectral diffu- 
sion, anh armo nicity effect ^^^ISOj^ phonon scattering from 
interface^^^^, and a modified phonon spectrunP3 Phe- 
nomenologically, we include radiative and pure dephasing 
mechanisms to broaden the ZPL, while radiative coupling 
to the cavity mode is treated in a self-consistent way. 

We seek to solve the Hamiltonian dynamics in the 
polaron fram^, where H' — > e^He~^, with P = 
a^a^ ^{b^q — bq). Closely following the notation of 

Wilson-Rae and Imamoglu^^-*, the polaron-transformed 
Hamiltonian can be written in terms of a modified sys- 
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which results in a shghtly different definition of the sys- 
tem Hamiltonian given in Eq. ( 2a I ; this modified system 



Hamiltonian is central to obtainin g result s in the nonper- 
turbative phonon-coupling regim^^HllIIIl, 

Figure 1 shows a schematic of an example QD-cavity 
system [Fig. 1(a)], as well as an energy level diagram 
showing the exciton levels [Fig. 1(b)], some of the phonon 
levels, and the cavity photon states. The effective 
exciton-cavity coupling term, (B) g, corresponds to a co- 
herent reduction of g due to acoustic phonon interactions 
[cf. Eq. pal)]. 



FIG. 1: (Color online) (a) Schematic of a semiconductor cav- 
ity used in cavity-QED (micropillar system), containing a cou- 
pled QD and driven by a coherent cw laser {r^x from the side). 
The fluorescence spectra of this QD-driven system is mea- 
sured via cavity emission (k) or QD emission (7). (b) Energy 
level diagram, where je) denotes the excited QD state, \g) 
denotes the ground QD state, and jO, 1, ■ • ■) show the cavity 
photon states. Also shown is the phonon reservoir (e.g., a 
bath of harmonic oscillator states treated at the level of the 
IBM) as green lines which is coupled to the QD excited state. 



tem, bath, and interaction part: 

H'g = h{A^L - Ap)cr+cr" + hAcLa^a + {B)Xg, (2a) 
H's=Y,f^^,b%, (2b) 



(2c) 



where B± are the coherent displacement operators of the 
phonon modes. 



B± = exp 



A, 



(3) 



and the fluctuation operators are defined through, Q — 
i(B+ + B_ - 2{B)) and C« = " B.). The 

thermally-averaged phonon displacement operators obey 
the relation'^ (B) = (B.) = (B), where 



(B) = exp 



dLo-^coi\i{phuj/2) 



(4) 



and J{oS) — A^(5(w — ujq) is the phonon spectral func- 
tion that will be explicitly defined later (in a more practi- 
cal continuum form). The polaron frequency shift, given 
by Ap = J2q^q/^q — /o dLdJ{uj)/uj, accounts for the 
renormalization of the QD exciton resonant frequency 
due to coupling to the phonon modes; however, we will 
assume that the polaron shift is implicitly included in 
ujx- The phonon-bath- modified system operators, Xg and 
Xu, are given by 



Xg = hg{a)(T + cr+a) + h7jx{(j + (t+). 



X,, 



ihgi^a^ a — a) a ) + ihrix{a^ 



(5a) 
(5b) 



B. Master equation 

The usual quantum optics approach to describe 
system-reservoir interactions is to derive a ME in a form 
suitable for numerical (or analytical) solution that can 
include photon interactions to any arbitrary order. The 
ME is usually derived using projection operator tech- 
niques such as the NZ approach or the TC approacfP^. 
The former technique leads to an integro-differential ME 
convoluted with the memory kernel which encapsulates 
the effects of the reservoir. However, the latter technique 
(TC) results in a ME which is local in time with time- 
dependent superoperators acting on the reduced density 
operatoiEsHIir 

The derivation of the ME for the reduced density oper- 
ator of the system in the NZ form starts with an equation 
that is written to second-order in the system-reservoir 
couplings (Born approximation in the interaction pic- 
ture. Defining the interaction Hamiltonian, 



i{H'g+H'^)t 



(6) 



then we obtain the following equation of motion for the 
reduced density operator. 



dm 



dt 



NZ 



1 



di'trB{[Hi(t), [Hi(t'),p(t')pB]]}, 

(7) 



where the tildes denote the operators in the interaction 
picture and 'tr' is the trace. A Markovian form of Eq. ([7|, 
to second order in the interaction terms, can be derived 
by approximating p(t') with p{t), so that 



dm 

dt 



TC 



1 



^ dt'trB{[Hi(t), [Hi(t'),p(t)/9B]]} 



(8) 



This TC form eliminates the p{t') dependence of the inte- 
gral on the right-hand side of the density operator equa- 
tion at all previous times'*^ However, some memory 
(time history) effects can still exist depending upon the 
form of the bath interaction. A further simplification to 
Eq. ^ is obtained by extending the upper limit of the 
integral to 00. This is a second-Markov approximation 
that yields a ME in the TC Markov form. In fact, this 
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latter (simpler) ME form can be further manipulated — 
under certain approximations — to produce an approxi- 
mate phonon ME of the standard Lindblad fornP3 we 
have previously demonstrated such an approach to study 
intensity power broadening. 

Substituting the interaction Hamiltonian [Eq. [6] into 
Eq. (|8|, gives 



Gmi~T) = ^^{t), and simplify Eq. Q as follows, 



dm 



dt 



TC 



--.2 



ra—g^u 



dt'iTB{(:rn{t)pBQm{t')}Xm{t)~p{t)X^{t') 



dt'tlB{Qm{t')pBC-m{t)}X^{t')p{t)X^m{t) 



+ / dt'iYB{pBU{t')Qm{t)}p{t)Xra{t')Xra{t) 

Jo 

(9) 



where Cm(t) = eJ-^B^CmG"^^^^ ■ The separation of the 
trace over the phonon variables is possible because we 
assume an initially separable density operator consisting 
of the system (reduced) density operator and the phonon 
bath density operator pB- We have also assumed that the 
full density operator [i.e, pBp{i)\ remains separable at all 
later times. We denote trB{ps^m(t)Cm(i')} = Gm{t — t') 
and tvB{pBC,m{t')C,m{t)} = Grn(t ' - t) , whcrc Gm{t) are 
the polaron Green function^^^I^ 



Gg{t)^ {Bf{cosh[cl>{t)]-l), (10a) 
GJt) = {Bf smhimi (10b) 



and the phonon correlation function (j){t) is defined as 



f°° J(uj) 

(j){t) = / duj-^ [coth(/3;iu;/2) cos(a;i) - i sin(w<)] . 
Jo ^ 

(11) 



The relevant phonon spectral function J(w) describes 
LA (longitudinal acoustic) phonon coupling via a 
deformation potential^, and is given by J(w) = 
apUJ^ exp(— w^/2w^) where ap is the exciton-phonon cou- 
pling strength and Wf, is the phonon cutoff frequency. 

Transforming back to the Schrodinger picture and 
changing variables from t — t' = r, we note that 



dp(t) 



dt 



TC 



= -,[H's.Pit)] 



]^ I dT J2 G^{T)X„,e-'"'s^"'X„,e'"'s-/f^p{t) 
^ /*^^ E G*^{r)X^p{t)e-^"'s-/^X„,e^"'s-/^ 

^ m—g,u 

\^ f dr G^{T)e-'"'s-"'X^e''''s-"'p{t)X^ 
-^/^^ E Gl,{T)p{t)e-^"'s-/^Xr^e^"s^/''^X^. 

^ m—g.u 

(12) 

This form of the TC ME is local in time and can be 
rewritten in compact form, 

m—g,u 

(Gm{r)[X„,,e-'''s-/'''X^e'"s-^'^p{t)]+U.c.) , 

(13) 

where we have added the C{p) superoperators to take into 
account additional dissipation processes that are neces- 
sary to obtain a complete description of the system-bath 
dynamic^^. We model these additional dissipation pro- 
cesses through the standard Lindblad operators: 

^{p) — -^{'ia'^ pcr^ — (T^cr^ p — p(T^a^) 
+ K,{2apa) — a) ap — pa) a) 

i 

+ — (2criipcrii - CTiicriip- pCTiiCTii), (14) 

where a\\ = ^ 2k is the cavity decay rate, 7' is the 

pure dephasing rate, and 72, = ^x{B)'^ is the radiative 
decay rate^^I^. Note that we include additional pure 
dephasing effects beyond the IBM, e.g., perhaps caused 
by quadratic QD-acoustic phonon and phonon-phonon 
interactions^^ with a rate 7'; this pure dephasing con- 
tribution to the ZPL broadening is al so kno wn to increase 
linearly with increasing temperature^^Mll 

The timescales associated with the phonon and radia- 
tive processes are substantially different Phonon pro- 
cesses happen at picosecond timescales whereas radiative 
timescales are smaller by (typically) at least an order of 
magnitudes). This results in a substantial simplification 
of the TC ME derived above. Extending the upper limit 
of the integral {t — ^ 00), Eq. (13) can be rewritten as 



dp{t) 



dt 



TC 



1 1 1'°° 

-,[H's,Pm+Cip)-^l dr E 



m—g,u 



{G^{T)[X^,e-^"s^/'''X^e'"s^/^p{t)]+R.c. 



(15) 
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This ME equation is now in a TC Markov form. The 
vahdity of extending the time hmits of the integration to 
c» depends on the ensuing system dynamics introduced 
by the external drive, and the exciton-cavity coupling 
(through g), with respect to the phonon bath relaxation 
time. For instance, a cavity-QED system driven by a 
short pulse may require a fully non-Markovian treatment 
to obtain a correct description of nonradiative processes. 
Non-Markovian effects within the polaron frame ME are 
likely to affect the population dynamics only at very 
short times and the difference between non-Markovian 
and Markovian dynamics vanish at longer timescale^^H, 
For the present system of interest (with cw drives), the 
Markov approximation is found to be an excellent ap- 
proximation as will be shown in the next section. 

Using similar steps as above, Eq. ([t]) leads to a ME of 
the NZ form that involves integ rating the system dynam- 
ics over all previous 

timePEl so that 

i^'^ m=g,u 

(16) 

This ME is non-local in time and generates equations of 
motion which are convoluted in time. 

We note that the validity of the Born approximation, 
used to derive the above MEs, depends on the physical 
parameters of the cavity-QED system. A general condi- 
tion for the validity of the second-order Born approxima- 
tion has been derived by McCutcheon and Nazij^S and is 

2 

given by %(1 — (B)'^) ^ 1. A similar condition for the 
validity of the Born approximation for the cavity-QED 
system can easily be derived with relevant numbers well 
below unity (since rix <C Wf, and g <^loi,, with = 1 meV 
for this work); this suggests that a second-order Born ap- 
proximation is perfectly valid for our present system. The 
main reason for this observation is that nonperturbative 
coupling effects beyond second-order Born are already 
captured within the polaron transform itself. 

One can also derive a ME describing the cavity-QED 
dynamics without using the polaron transformatiorl^. 
One common approach is to use a standard weak-coupling 
theory which only includes exciton-phonon interactions to 
second-order and is known to breakdown for higher tem- 
peratures and strong exciton-phonon coupling^^Il with 
cavity coupling, thes e weak phonon-coupling failures can 
become even stronger'^®'^^'. A more general approach us- 
ing a variationally-optimized ME has also been described 
by McCutcheon et al^^. The variational ME includes 
nonperturbative phonon effects and multiphonon pro- 
cesses over a wide range of parameters (especially im- 
portant for very large driving fields). Glassl et al.^'' 
have gone beyond a polaron picture as well, and car- 
ried out numerically-exact real-time path integrals and 
a fourth-order correlation expansion to model the long- 
time dynamics and the stationary nonequilibrium state 



of an optically driven QD (coupled to acoustic phonons) . 
For the purposes of our paper, as noted above, the po- 
laron transformation is rigorously valid as the coherent 
laser pump rates and the exciton-cavity coupling rate are 
much smaller than the phonon cutoff frequency ^^'l More- 
over, a clear advantage of our approach is the ease with 
which cavity coupling is included, and, for example, the 
well known Jaynes-Cummings model is formally recov- 
ered with phonon coupling switched ofiP^. 

With coherent pumping, we define the incoherent flu- 
orescence spectra associated with cavity mode emission 
(S'c), and QD emission {Sx^^, as follows: 



Sc{oj) cc — lim Re 



dr {{a\t + T)a{t)) 



{a{t)){aHt)))e 



Sxi^jj) oc — lim Re 

7f f oo 



(17) 



dT{{<j+{t + T)a-{t)) 



((T+(<))(a-(t)))e"*'(^)e^("^-'^)^ 



(18) 



We highlight the phase term e"^ ^■'^^ appearing in the def- 
inition of the exciton spectra, S^iuj), which arises due 
to the polaron transformation and accounts for the de- 
phasing of the exciton coherence in the phonon reservoir. 
In the polaron frame, the exciton autocorrelation func- 
tion {a^{t + T)a^{t)) is transformed to {a'^{t + T)B^{t + 
t)(j~ {t)B-{t)) . The trace over the phonon degrees of 
freedom and the QD Hilbert space is easily evaluated as 
the Hilbert space is separable in the polaron frame. 

It is also useful to define the QD susceptibility (or 
polarizability), x(w), that describes the linear response 
of the QD in the presence of phonon coupling and phe- 
nomenological dissipation only (i.e., not coupled to the 
cavity mode, g — 0, and driven by a linear optical field). 
This bare polarizability function is obtained from 



X(w) =x'(w)+«x"(w) 



X(c^) cx / dT{n-{T)a+{Q))e 
Jo 



(19) 



where we now use a rotating frame with respect to uj^. 
This simple looking form for x(a;), with the phonon con- 
tribution described by the phonon correlation function 
(f>{t), is due to the polaron transformation. To obtain 
the above QD response function, the expectation is eval- 
uated in the absence of any exciton-cavity coupling, and 
we excite the QD exciton with a linear (weak) inco- 
herent pump process, via the Lindblad superoperator: 
Px{2a^ pa~ —a~a^p — pa~a^) with a suitably small P^. 
The linear absorption spectra is defined through x" {^)- 
It is interesting to recognize that the phonon correla- 
tion function 4>{t) is explicitly present in the definition of 
the Sx and x('^)- The dynamics of equations we solve is 
effectively Markovian in the polaron frame, which is much 
easier to work with from a numerical viewpoint. How- 
ever, the transformation back to the original frame rein- 
troduces the non-Markovian phonon dynamics through 
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the phonon correlation function (j){t). Note also that 
the cavity spectrum does not depend directly on (/'(i), 
so the cavity spectrum may not depend at all on any 
non-Makovian dynamics (at least from a mathematical 
viewpoint, and if one works in a polaron frame). 

The fluorescence spectra are calculated by the numer- 
ical solution of the ME driven by a steady-state pump 
r]x. The ME is solved in a Jaynes-Cummings basis 
with states |0), \1L), \1U), \2L), \2U) and so on. De- 
noting the photon states \n), with n = 0,1,2,..., and 
exciton states \e/g), the Jaynes-Cummings basis states 
are given by |0) = \g) \0),\1L) - 1/V2{\g) |1) - |e) |0)), 
\1U) = l/V^Qg) |1) + |e) |0)), etc. The system is ini- 
tialized with the exciton and cavity in the ground state 
which then evolve, in the presence of a cw pump, to a 
dynamical steady state. Using the quantum regression 
formulcP^, we subsequently compute the relevant two- 
time correlation function whose decay is initialized by 
the steady state population and coherences. The calcu- 
lations are performed, in part, using the quantum op- 
tics toolbox by TarfSSl. We note that, in general, it is 
vital to include multiphoton effects to get the correct 
(numerically converged in the photon number) fluores- 
cence spectra; though a one photon basis is (obviously) 
enough to obtain the linear absorption. We have verified 
the need for multiphoton effects elsewhere and find that 
the truncation of the hierarchy of the photon states to 
two-photon-correlations (two photon states) is sufficient 
(and necessary) for exciton driving^. A broader discus- 
sion of the role of multiphoton processes, especially in 
the context of intensity power broadening, is presented 
in Ref. ESI 



III. NUMERICAL RESULTS 

In this section, we present and discuss our numer- 
ical results. For all calculations that follow, we use 
the following phonon parameters: w;, = 1 meV and 
ap/(27r)^ = 0. 06 ps^, which are typical numbers for 
InAs/GaAs QDJ^HEUmI^ Pqj. ^j^g cavity decay rate, we 
choose K = 50 /xeV, and for the exciton decays, — 
1 /ieV and 7' = 2 /ieV. For simplicity, we will fix the 
pure dephasing rate below, but thi s rate is also known 
to increase linearly with temperatur^^HlSl -pj^^g other pa- 
rameters are generally varied or given in the figure cap- 
tions. 



A. Linear absorption spectra 

In Fig. [2] we show the linear absorption spectra of the 
QD exciton, x"(w), for two different temperatures of the 
phonon bath. The absorption spectra in the vicinity of 
the QD resonance is dominated by the IBM spectral func- 
tion which allows the dot to absorb photons by phonon- 
assisted processes. In the absence of any phonon coupling 
(blue dashed lines in Fig.[2|, the absorption spectrum is 
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FIG. 2: (Color online) Linear absorption spectra of the 
QD, x"{^)j for two different phonon bath temperatures: (a) 
T = 4 K and (b) T = 10 K. The red (soUd) line includes 
phonon coupling and blue (dashed) hne has no phonon cou- 
pling. The QD decay parameters jx ~ 1 /^eV and j' — 2 fieV 
are used in both plots, which give the broadening of the ZPL. 
The asymmetry in the absorption spectra due to the phonon 
reservoir is clearly seen, especially at lower temperatures. 



a simple Lorentzian whose width is determined by radia- 
tive broadening and pure dephasing (the ZPL linewidth) . 
The asymmetry in the absorption spectra in the presence 
of phonons is clearly seen at low temperatures, which 
is caused by different probabilities to absorb or emit a 
phonon; i.e., the latter process is more likely at lower 
temperatures. 



B. TC vs NZ master equation solutions: 
population dynamics and Mollow triplet spectra 

In Fig. |3] we compare the population dynamics gener- 
ated by the NZ ME [Eq. the TC ME [Eq. ([l3|] 
and the TC Markov ME [Eq. ([l5])]. We excite res- 
onantly (i.e., ujl = Wx), with the cavity detuned to 
lower energies (a;^, — w^, = 0.28 meV), using the exci- 
ton pump rate, 7]^ = 0.05 meV. As can be seen from 



the plots, all three ME forms produce identical results 
for the evolution of the cavity and QD populations, de- 
fined as ndt) = {a^a{t)) and n^(t) = {a^a^{t)). We also 
compare the cavity-emitted fluorescence spectra (Sc) as 
shown in Fig. |4] (defined in the previous section) , gener- 
ated by the respective MEs, which are again found to be 
identical (or, at least, indistinguishable). Hence, in the 
remainder of the paper we use the TC Markov form of 
the ME as it is much easier to work with and is justified 
here to be an excellent approximation. Note that for the 
bath temperature of 10 K, (B) = 0.85, which reduces 
the effect of both rjc and g, and causes a smaller triplet 
splitting. 

For numerical convenience, the calculations in this sub- 




200 400 600 800 1000 

Time (ps) 

FIG. 3; (Color online) (a) Comparison of the population 
dynamics generated by the NZ and the TC form of the 
ME, both to second-order in the system-phonon-resorvoir 
coupling. The upper blue (dark solid) line is the QD exciton 
population nx{t) using the NZ solution, and the orange (light 
dashed) line is the TC ME solution. We also plot the cavity 
mode population nc(t). The grey (light solid) lines show the 
exciton and the cavity mode populations in the absence of 
any microscopic phonon coupling (i.e., only including ZPL 
broadening and cavity coupling). In (b) we compare the 
population dynamics of the TC ME and the TC Markov ME. 
The blue (solid) line shows the exciton population iix (t) using 
the TC ME, while the orange (dashed) line is the dynamics 
generated by the TC Markov ME form. All the plots were 
obtained by truncating the hierarchy of equations of motion 
at the one-photon limit, and using T = 10 K, 7^; = 1 /ieV, 
K = 50 neV, 7' = 2 yueV, 5 = 20 /xeV, and rj^ = 0.05 meV. 
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FIG. 4: (Color online) (a) Corresponding cavity-emitted res- 
onance fluorescence spectra from Fig. 3. The blue (solid) line 
is the cavity-emitted fluorescence spectrum (Sc) computed us- 
ing the NZ form and the orange (dashed) line is the spectrum 
generated by the TC ME. The grey solid lines show the spec- 
trum in the absence of any phonon coupling (apart from the 
ZPL) . In (b) we compare the cavity-emitted spectra generated 
by the TC and TC Markov ME resuhs. 



mode as the cavity-QD coupling is further strengthened 
by phonon emission processes. At a detuning of ± 1 meV, 
the red-detuned cavity peak is found to be about dou- 
ble that of the blue detuned cavity mode peak, which is 



section (shown in Fig. [s] and Fig. |4]) were obtained by 
truncating the hierarchy of equations of motion at the 
one photon limit because of the computational complex- 
ity of solving the NZ ME form. However, all the results 
discussed in the subsequent sections use the TC Markov 
ME with the required (i.e., numerically converged) num- 
ber of multiphoton processes; so the fluorescence spec- 
trum will be seen to qualitatively change, especially near 
the cavity resonance (see Ref . for more details). 



C. QD-cavity detuning dependence of the 
phonon-dressed Mollow triplet 

In Fig. [5] we study the detuning dependence of Sc for 
different QD-cavity detunings. The fluorescence spec- 
tra are clearly asymmetric with respect to the cavity 
mode detuning; this behavior is drastically different to 
simple Lorentzian-decay models which would yield sym- 
metric Mollow spectra, i.e., with respect to the exciton- 
cavity detuning. Cavity emission is considerably en- 
hanced when the QD is detuned to the right of the cavity 
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FIG. 5: (Color online) Detuning dependence of Sc for different 
exciton-cavity detunings. The other parameters are the same 
as in Fig. |3j again with the phonon bath temperature at T = 
10 K. In (a) the red lines denote Sc when the cavity is detuned 
to the left of the QD (exciton). In (b) the blue lines denote 
the same spectra when the cavity is detuned to the right of 
the QD. Note that all the spectra are normalized relative to 
the first spectrum on the top left. 
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consistent with experimental data using incoherent ex- 
citation, e.g., as shown by Hennessy et al^. Also note 
that the cavity-mode peak initially increases, relative to 
the MoUow triplet, as the QD exciton is red detuned from 
the exciton; but eventually decreases for sufficiently large 
QD-cavity detunings. This suggests that the coupling 
between the QD and the cavity is optimum at a par- 
ticular detuning and is determined by the QD-phonon 
coupling strength and the intrinsic spectral properties of 
the phonon reservoir. This observation is significant for 
experimentalists as it suggests that off-resonance cavity- 
QED physics is best studied at a particular QD-cavity 
detuning where the mode emission is maximized. Note 
that the emission at the cavity mode frequency is even- 
tually suppressed at very large QD detunings A^^: — ±5 
meV, which is in agreement with recent experimental ob- 
servations on site-controlled QDs in cavitieJ^ however, 
these measurements were again carried out using inco- 
herent excitation, so not in the domain of resonance flu- 
orescence. The maximum coupling between the QD and 
the cavity also depends on temperature, e.g., the QD- 
cavity detuning at which the QD-cavity mode coupling 
is maximum increases with temperature. 

It is informative to note that the detuning dependence 
of the spectra closely follows the spectral characteris- 
tics of the phonon correlation function. An intuitive 
way to understand this detuning dependence is using 
effective phonon rates in a Lindblad-type descriptiorP^, 
which are similar to the cavity-feeding rates introduced 
in Refs. 39 40, These phonon scattering rates are related 
to the Fourier component of the phonon correlation func- 
tion at the QD-cavity mode detuning. The same type of 
feature can be seen in the asymmetry of the sidebands 
of the vacuum-Rabi doublet in an incoherently pumped 
strongly coupled QD-cavity systerrP^l the asymmetry of 
the doublet depends closely on the spectral character- 
istics of the phonon correlation function. Hence, these 
features in the fluorescence spectra and linear absorption 
reveal subtle details of the underlying phonon reservoir, 
which can be probed as a function of the exciton-cavity 
detuning. The advantage of coherently excited systems 
is that they are much cleaner to study, since for inco- 
herent excitation one never really knows the details of 
the pumping scenario (which may involve a number of 
exciton states and contributions from different QDs). 



D. Temperature dependence of the resonance 
fluorescence spectra 



frequency increases with temperature. In Fig. |6jb), we 
plot the QD exciton spectra, where the fuU-polaron QD 
spectra is now plotted as the blue (dark) lines. 

The one-phonon limit of the full-polaron ME is ob- 
tained by a perturbative expansion of the relevant 
phonon correlation functions to lowest order; namely, 
we approximate Ggij) ~ and Gu{t) — (t>{t). The 
TC Markov ME is then solved numerically for the above 
choice of phonon correlation functions to study the flu- 
orescence spectra in the one-phonon limit. In the weak 
phonon coupling theory, (B) ~ 1, suggesting that one- 
phonon processes do not renormalize the Rabi frequency 
of the cw-laser drive, rj^, or the QD-cavity coupling con- 
stant g. The discrepancy between the one-phonon and 
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In Fig. [6] we investigate the temperature dependence of 
the resonance fluorescence spectra and also analyse the 
role of one-phonon scattering versus multiphonon scatter- 
ing. In Fig. |6|^a), we show the red-detuned cavity mode, 
where the grey (solid) lines denote the one-phonon limit 
of the full polaron cavity-emitted spectra which is shown 
with the red (dark) line. We observe that the emission 
at the cavity mode relative to the emission at the QD 



FIG. 6: (Color online) Temperature dependence of the (a) 
cavity-emitted resonance fluorescence spectra (red, dark line) 
and (b) exciton-emitted resonance fluorescence spectra (blue, 
dark line) for different temperatures of the phonon reservoir. 
Also shown is the one-phonon calculations (grey, light curves). 
The QD-cavity detuning, Acx = —0.5 meV, with the phonon 
reservoh temperatures: T = 4 K, T = 10 K, T = 20 K, and 
T = 40 K. The other parameters are the same as in Fig. [3] 
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FIG. 7: (Color online) Resonance fluores- 
cence spectra for various cw-laser drives 
and two different exciton-cavity coupling 
g. The phonon reservoir at T = 10 K, and 
we use Acx = —0.5 meV, ~ I f^eV, 
K = 50 /xeV, and — 2 /leV. In Figs, (a) 
and (b) we plot Sc and Sx for exciton- 
cavity coupling g = 20 fieV with various 
cw-laser driving rjx which are shown as red 
(dark) lines and blue (dark) lines, respec- 
tively (left and right panels) . Also plotted 
are the grey (solid) lines which are the cor- 
responding fluorescence spectra in the ab- 
sence of any phonon coupling (apart from 
the inclusion of the ZPL broadening) . In 
Figs, (c) and (d) we plot Sc and Sx with 
the larger g = 50 fieV using similar pa- 
rameters as in Figs, (a) and (b) . 



full-polaron solution of the TC ME increases with tem- 
perature, and the one-phonon solution is found to over- 
estimate the EID in comparison with the full polaron 
solution. These results are consistent with the findings 
of McCutcheon and NaziJ^S!^ who adopt a polaron ME for 
analysing QD Rabi oscillations in the presence of acoustic 
phonon scattering (with no cavity). Note that since the 
onc-phonon truncation for the full polaron theory does 
not renormalize the laser Rabi frequency, the locations 
of the triplets are identical to that for the case of no 
phonons. Hence a full polaron theory can be essential to 
find the correct center frequencies of the MoUow triplet 
resonances. Moreover, the discrepancy between the side- 
band locations increases with temperature. However, the 
weak phonon theory can obtain accurate simulations of 
the cavity resonance, depending upon the exciton-cavity 
detuning. Without any acoustic phonon coupling, the 
Mollow spectra are significantly different, e.g., there is no 
EID and there is negligible oscillator strength at the cav- 
ity mode. These results suggest that a simple atomiclike 
ME would drastically fail in the present semiconductor 
cavity-QED excitation regime. 



E. Resonance fluorescence spectra versus exciton 
pump rate 

In Fig. [T] we display the normalized resonance fluores- 
cence spectra for various cw laser drives and dot-cavity 
couplings, using a fixed QD-cavity detuning. The tem- 
perature of the bath is at T = 10 K. From Fig. and 
Fig.[7|^c), we observe that the relative emission at the cav- 
ity mode frequency increases with increasing drive which 
is caused by the enhanced coupling between the cavity 



mode and the QD, as the low energy triplet component 
is pushed closer to the cavity mode frequency. The Mol- 
low triplet peak which is proximal to the cavity mode 
has a higher spectral weight than the distant component 
which is due to stronger relative coupling with the cav- 
ity mode. However, we note that the asymmetry in the 
spectral weights decreases for larger g as both the proxi- 
mal and the distant component of the triplet are more or 
less equally coupled to the cavity mode. For increasing 
pump strengths, it is clear that the linewidth broadening 
is larger in the presence of phonons. Further increased 
broadening of the Mollow sideband occurs with increas- 
ing 

Apart from the lack of any emission at the mode fre- 
quency in the QD emission spectra, the asymmetry in the 
Mollow triplets in Sx is substantially smaller than Sc- 
Furthermore, the high energy component of the triplet 
has a larger relative spectral weight than the low energy 
component. This asymmetry in the sidebands spectral 
weights increase with increasing g and cw-laser drive rjx ■ 
Also, the asymmetry of the sidebands are not mainly due 
to phonon-induced effects as they are present even in the 
absence of phonons; thus these features are primarily cav- 
ity coupling effects and are easy to understand in terms of 
relative coupling with the cavity mode. The low-energy 
component of the triplet is closer to the cavity mode and 
couples more strongly with it. 



F. Fluorescence spectra with an off-resonance cw 
drive 

In Fig. [8] we study Sc and S^ for different exciton-laser 
detunings, using two different exciton-cavity couplings, g. 
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FIG. 8: (Color online) Similar to Fig. [t] 
but with rix fixed at 50 /leV and for dn- 
ferent exciton-cw laser detunings Alx- In 
Figs, (a) and (b) we plot the resonance 
fluorescence spectra (5*^ and Sx, respec- 
tively) for QD-cavity coupling g — 20 fieV 
with various exciton-laser detunings re- 
spectively plotted as red (solid) lines and 
blue (solid) lines. Also shown are the flu- 
orescence spectra in the absence of any 
phonon coupling as grey-solid lines (apart 
from the ZPL). In Figs, (c) and (d) we plot 
the fluorescence spectra for exciton-cavity 
coupling g = 50 /xeV with identical detun- 
ings and cw laser drives as in Figs, (a) and 
(b). 



Again we fix the phonon reservoir at T = 10 K. We first 
note that, unlike the on-resonance spectra (see Fig. [?]), 
the relative spectral weights of the MoUow triplets are 
also now determined by A/^j. (the exciton-cw laser de- 
tuning). The high-energy component of the MoUow 
triplet has larger relative spectral weight if the cw-laser 
is detuned to the left of the exciton (lower energies). 
Conversely, the low-energy component of the triplet has 
larger relative spectral weight if the cw-laser is detuned 
to the right of the exciton (higher energies). This ob- 
servation is also consisent with the features of Sx- An 
ofF-resonant cw-laser drive shifts the entire triplet with a 
corresponding shift in the location of the sidebands. A 
positive exciton-cw laser detuning shifts the triplet struc- 
ture further away from the cavity mode but increases the 
spectral weight in the low energy sideband as it is pushed 
closer to the bare exciton resonance. A negative exciton- 
cw laser detuning shifts the triplet structure towards the 
cavity mode, and increases the spectral weight in the high 
energy sideband as it is pushed closer to the bare exciton 
resonance. Also note the asymmetry in the low-energy 
sideband of the triplet in the top panel and the high- 
energy sideband in the bottom panel in both Fig. [sf^a) 
and Fig. [sjjc) which is explained by the difference in the 
interaction with the cavity owing to different effective 
spectral detuning between the respective sideband and 
the cavity mode. The asymmetry in the spectral weights 
is suppressed with increasing QD-cavity coupling as seen 
by comparing, e.g., the relative weights in Fig. [sj^a) and 
Fig. [sjjc) . The enhanced exciton-cavity interaction over- 
rides, to some extent, the effects of the spectral shift 
induced by the off-resonant drive. Finally, we note that 
phonon coupling enhances the asymmetry in the relative 
spectral weights of the MoUow triplets due to enhanced 



exciton coupling and reduced sideband splitting. 



G. FWHM (linewidth) of the Mollow sidebands 

We next study the broadening of the Mollow sidebands 
for various excitation regimes. 

In Fig. |9] we show the numerically extracted FWHM 
(full width at half maximum) of the Mollow triplet peaks 
of the cavity-emitted spectrum (Sc) for various exciton- 
cavity detunings and exciton-cavity couplings as a func- 
tion of the square of the Rabi frequency rj^- Clear sig- 
natures of phonon-induced EID are demonstrated by the 
dependence of the FWHM on the square of the Rabi fre- 
quency; for smaller exciton-cavity couplings [Fig. [9f^a) 
and Fig. [9f^d)], this dependence is approximately lin- 
ear (for the chosen exciton-cavity detuning). These two 
plots have opposite QD-cavity detunings, and EID is 
marginally larger when the cavity mode is detuned to 
the left of the QD due to increased phonon-induced cou- 
pling between the exciton and the cavity mode. This 
general trend of EID increasing linearly with the square 
of the Rabi frequency is in good agreement with recent 
experiments performed in semiconductor micropillarJl^l^ 
where ^oxp ~ 16 /ieV [cf. Fig. [9][d)]. Increasing the 
QD-cavity coupling g also increases the FWHM of the 
sidebands due to stronger coupling between the QD and 
the cavity and enhanced emission into the cavity mode 
[Fig.[9];b) and Fig.[9];e)]. The dependence at higher QD- 
cavity couplings is nonlinear especially at low drives [see 
Figs.[9]jc) and Fig. [ojf)]; we also find a linear increase in 
the FWHM as a function of the square of cw laser drive 
even in the absence of phonon coupling, primarily for 
large QD-cavity couplings. This effect can be attributed 
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FIG. 9: (Color online) Numerically extracted 
FWHM of the on-resonance MoUow triplet spec- 
trum emitted via the cavity mode {Sc) for various 
exciton-cavity detunings and g, as a function of the 
square of the Rabi frequency rjx- In Figs, (a)-(c) 
we plot the FWHM of the sidebands of the fluo- 
rescence spectrum for Acx = 0.3 meV. The orange 
(light) data is with phonon coupling and the blue 
(dark) data is for no phonon coupling. The crosses 
denote the high frequency and the circles denote 
the low frequency component of the MoUow triplet. 
Also plotted as a solid line with corresponding color 
is the best linear fit to the data. In Figs, (d)-(f) 
we show the FWHM of the sidebands of the fluores- 
cence spectrum for Acx = —0.3 meV, using identical 
parameters as in the above. 
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FIG. 10: (Color online) Numerically extracted 
FWHM of the off-resonance Mollow triplet spec- 
trum emitted via the cavity mode (Sc) for different 
exciton-cavity detunings and g, as a function of the 
QD-laser detuning A^l- In Figs, (a)-(c) we plot the 
FWHM of the sidebands of the fluorescence spec- 
trum in the absence of phonons with the exciton- 
cavity couplings: g = 20 fieV, g = 50 fieV and 
g = 80 fieV. The blue data is with Acx = 0.3 meV 
and the red data is for Acx ~ —0.3 meV. The crosses 
denote the high frequency and the circles denote 
the low frequency component of the Mollow triplet. 
The triangles are the average of the two sideband 
FWHM. In Figs, (d)-(f) we study the FWHM of 
the sidebands of the fluorescence spectrum, with 
phonon scattering at T = 10 K. 



to enhanced emission into the cavity mode which is en- 
hanced with increasing drive and broadens the triplet. 
We also observe a notable difference in the broadenings of 
the low-energy and high-energy components of the triplet 
in the presence of phonons. The triplet component prox- 
imal to the cavity mode shows larger broadening which 
increases with g. 

In Fig. [To] we study the numerically extracted FWHM 
of the ofF-resonance Mollow triplet as a function of the 
exciton-laser detuning A^x- The coupling to electron- 
phonon scattering is seen to introduce additional broad- 
ening of the sidebands, as expected, due to FID; this 
can be seen from the differences in magnitudes between 
the plots in the top panel and the bottom panel. Also 
note that the two sidebands of the triplet are broadened 
by different amounts which is governed primarily by the 
proximity of the sideband to the cavity mode. In the 



absence of any phonon coupling, the Mollow component 
which is closer to the cavity mode is broadened more 
due to enhanced emission into the cavity. Moreover, the 
difference in the broadening of the high energy and the 
low energy components of the triplet increases with tem- 
perature. We also observe a noticeable difference in the 
broadening of the sidebands in the presence of phonons. 
We find an enhancement in the FWHM broadening when 
the cw laser drive is spectrally located between the cavity 
mode and the exciton mode. In contrast, there is sub- 
stantial suppression of FWHM if the cw laser is detuned 
in the other direction relative to the exciton-cavity detun- 
ing. This suggest that phonons enhance the broadening 
process in the first scenario as it allows for an additional 
mechanism for enhanced emission into the cavity mode. 

Recent experiments in semiconductor micropillars'iSl 
have noted a distinct spectral narrowing of both the 
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high energy and the low energy sideband of the MoUow 
triplet with increasing exciton-cw laser detuning. This 
behaviour is quite contrary to what is obtained from a 
phonon-based EID model. However, a decrease in the 
Rabi frequency was also noticed, experimentally, as a 
function of detuninjl^l jj^ addition, more recent experi- 
mental data obtained for weakly coupled QD-cavity sys- 
tems find the opposite dependence on the exciton-laser 
detuning wherein the FWHM of the off-resonance Mol- 
low triplet is found to increase with the QD-laser detun- 
in^SH, Part of this anomaly can perhaps be attributed 
to additional detuning-dependent pure dephasing mech- 
anism^i^. Further investigations and a closer collabora- 
tion between experiments and theory are likely required 
to understand these conflicting trends in more detail. 

IV. CONCLUSIONS 

We have presented a detailed theoretical analysis of 
the fluorescence spectra from a coherently driven cavity- 
QED system. Fundamentally distinctive features of the 
semiconductor QD system are found because of acoustic 
phonon scattering in the solid state environment, which 
we have modelled at a microscopic level. The nonpertur- 
bative treatment of the effects of QD exciton, acoustic 
phonons, cavity and the cw laser field on the coherent 
part of the Hamiltonian, was accomplished by a suitable 
choice of the decomposition of the system Hamiltonian. 
A polaron transformation was adopted that retains co- 
herent phonon coupling to all orders, while introducing 
a suitable (modified) system-bath-interaction with which 
to perform perturbation theory that is accurate to sec- 
ond order. We derived a TC ME in the polaron frame to 
study the fluorescence spectra, and we showed that this 
ME produces identical numerical results to a NZ ME 
solution — the latter being significantly more difficult to 
solve numerically. We also showed that a Markovian ME 
is perfectly valid for the cw pumping scenarios investi- 
gated in this paper, if one uses a polaron frame. 

We then calculated the Mollow spectra as a func- 
tion of exciton-cavity detuning and found that electron- 
phonon-scattering causes EID and significant coupling 
between the exciton and cavity mode, throughout the 
entire phonon spectral function (easily spanning 5 meV) . 
At a phonon bath temperature of T = 10 K, the cav- 
ity resonance is found to go through a maximum at 
a specific detuning, when the cavity is red shifted by 
about 1.4 meV. The resonance spectra were found to be 
strongly asymmetric with respect to exciton-cavity de- 
tuning. Next, we showed the role of pump excitation 



(rjx) and exciton-cavity coupling rate {g), both of which 
are seen to have a significant impact due to phonon- 
induced scattering processes. Both the cavity-emitted 
spectra (Sc) and exciton-cmittcd spectra (Sx) were cal- 
culated and these were found to have quite different spec- 
tral forms, especially near the cavity mode resonance. 
We subsequently investigated the spectral FWHM of the 
Mollow triplet sidebands for various exciton-cavity de- 
tunings and exciton-cavity couplings. We analysed the 
Mollow triplet spectra as a function of the square of the 
Rabi frequency, and obtained unambiguous signatures of 
phonon-induced EID in accordance with recent experi- 
ments. We also studied the FWHM of the off-resonance 
Mollow triplet as a function of the exciton-laser detun- 
ing. The interaction of the QD with the phonon reservoir 
was found to introduce pronounced features in the fluo- 
rescence spectra exemplified by emission at the cavity 
mode and additional spectral broadening effects. With 
a systematic detuning between the driving laser and the 
exciton, our theoretical results predict further broaden- 
ing which is a prediction that is in disagreement with re- 
cent experiments in semiconductor micropillars^^ — which 
suggest a spectral narrowing of sidebands with increas- 
ing exciton-laser detuning. Nevertheless, as mentioned 
above, recent experiments also show a spectral broaden- 
ing with exciton-laser detunin^^. The possible role of 
detuning-dependent pure dephasing also merits further 
study in this context, as relative magnitudes of pure de- 
phasing and the QD radiative decay rate determine the 
curvature of the FWHM plot for the off-resonant fluo- 
rescence sideband as a function of QD-laser detuning. 
Additional contributions from phonon-damped processes 
and higher order QD-phonon interactions may further 
influence the underlying physics of the laser-detuning- 
dependent EID. 

We highlight that the semiconductor-specific features 
highlighted above are quite distinct to the physics of 
atomic cavity-QED and in general the application of 
an atomiclike MEs to these systems may drastically 
fail to produce the qualitative experimental features of 
semiconductor-based cavity-QED systems. Although we 
have specialized our studies to an exciton-driven system, 
it is straightforward to also study a cavity-driven system 
as well, e.g., see Ref. l55l 
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